Asteroseismology of Procyon A with SARG at TNG* 
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Abstract. We present high precision radial velocity measurements on the F5 IV star a CMi obtained by the SARG spectrograph 
at TNG (Telescopio Nazionale Galileo) exploiting the iodine cell technique. The time series of about 950 spectra of Procyon 
A taken during 6 observation nights are affected by an individual error of 1.3ms . Thanks to the iodine cell technique, the 
spectrogr aph contribution to the Doppler shift measurement error is quite negligible and our err or is do minated by the photon 
statistics iBrown et alll994) . An excess of power between 0.5 and 1.5 mHz, detected also bv lMartic et al. I l2004|) h as been 
found. We determined a large sep aration frequency Av p = 56+2 //Hz, consistent with both theoretical estimates I Chabover et al. 
1999) and previous observations (Ma rtic et al. l2004h . 
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1. Introduction 

Procyon A (a CMi, HR 2943, HD61421) for its proximity 
and bright ness has already attracted attention of stellar seis- 
mologists JGuenther & Demaraue Ill993t iBarban et al. l fl999: 
IChaboveretal.lll999HMarticetal.lll999l) . It is a F5 IV star 
(with V= 0.363) at a distance of 3.53 pc in a 40-year period 
visual binary system; the companion is a white dwarf more 
than 10 mag fainter than the primary. Adopting the very pre- 
cisely measured parallax by Hipparcos, II = 285.93+0.88 mas, 
Prieto et al. (2002) derived a mass of 1 .42 + 0.06 M , a radius 
■ of R/R Q = 2.071 + 0.02 and a gravity of log g = 3.96 + 0.02. 
These parameters cannot unambigously determine the evo- 
lutionary status of the star, which could be either in the 
core hydrogen-burning phase or in the more advanced hy- 
drogen shell-burning phase. Although in principle seismology 
would help to est ablish the evolutionary status of Procyon A 
jPi Mauro II2004I) . the nature of its pulsational spectrum is a 
subject of intense debate. 

In fact, the exc ess of power around 1 mHz found by 
Mai1ic_et_ajJi2004) seems to have a stellar origin and is con- 
sistent with a p-mode comb-like pattern with a large fre- 
quency separation of about 54 pHz. Howeve r, data from the 
recently launched Canadian MOST satellite ( Matt hews et al. I 
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2004) show no si gnificant power excess in the same spec- 
tral region. While Matthews et al. I (120041 suggested that the 
ground-based detection of p-modes may be an artefact caused 
by a combinat ion of stellar noise, and da ta sampling and analy- 
sis, [christensei^Dailgi^^ suggest that the 
most likely explanation for the null detection seems to be a 
dominating non-stellar noise source in the MOST data. 

The aim of this letter is to discuss the observations of 
Procyon A carried out using the high-resolution spectrograph 
SARG mounted at TNG (Telescopio Nazionale Galileo), aimed 
at detecting stellar oscillations by means of high-precision ra- 
dial velocity measurements. Our findings confirm the presence 
of an excess of power around the 1 mHz region, where we de- 
tect several frequencies whose distribution is consistent with 
a p-mode spectral pattern. In addition, our results show the 
high efficiency of the SARG spectrograph in detecting solar- 
like oscillations in stars, being an ideal instrument for multisite 
ground-based observing campaigns. 



2. Observations 

SARG is a high resolution cross dispersed echelle spectrograph 
(Grat ton et al. fc oOl) which offers both single object and long 
slit (up to 26 arcsec) observing modes covering a spectral range 




Fig. 1. A segment of the echelle spectrum of Procyon obtained 
by SARG. Top panel: iodine absorption lines are superimposed 
to the Procyon lines. Bottom panel: the Procyon spectrum with- 
out iodine absorption features. This spectrum has been used as 
the template spectrum in the determination of radial velocity 
from star plus iodine spectra 

Table 1. Observing log: date, number of spectra, signal-to- 
noise ratio and seeing. 



Date 


Sp. number 


SNR 


Seeing (as) 


2001/01/02 


160 


325 +4 


1.04 ±0.02 


2001/01/03 


149 


248 + 4 


0.81 +0.01 


2001/01/04 


181 


292 + 3 


0.82 + 0.01 


2001/01/06 


181 


295 + 5 


0.93 + 0.01 


2001/01/08 


125 


291 +5 


0.89 + 0.02 


2001/01/09 


153 


328 + 4 


0.87 ± 0.01 



Fig. 2. Radial velocities of Procyon obtained on 2001/01/06. 
The data are affected by an internal error of 1.48 ms -1 and an 
r.m.s. of 3.5 m s _1 . The 21-minute pulsation is also reported. 




Frequency (mHz) 



from A = 370 nm up to about lOOOnm, with a resolution rang- 
ing from R = 29, 000 up to R = 164, 000. 

Our spectra were obtained at R — 144, 000 in the wave- 
length range between 462 and 792 nm. The iodine cell covers 
only the blue part (see Fig. 1) of the spectrum (from 462 up 
to 620 nm) which has been used for measuring stellar Doppler 
shifts. The red part of the echelle spectrum of Procyon A has 
been used for measuring equivalent width of absorption lines 
sensitive to temperature. Here we present the results concern- 
ing the blue part of the spectrum while the analysis of the red 
one will be the topic of a future work. During the observations 
(the first SARG scientific run after commissioning and test- 
ing phases) we collected about 950 spectra with high signal- 
to-noise ratio and a mean exposure time of about 10 s. Due to 
weather and technical conditions (see Table 1 for more details) 
a few gaps are present in the time series. 

3. Data analysis 

Radial velocities have been obtained using the AUSTRAL code 
fcndletal.ll2OO0l) which models instrumental profile, star and 
iodine cell spectra in order to measure Doppler shifts. Since 
the instrument profile changes along a spectral order, the spec- 
trum was divided in segments (chunks) approximately 100 pix- 
els long (~ 2 A), each one modeled independently. About 400 



Fig. 3. The Scargle-Lomb periodogram of the data. An excess 
of power around 1 mHz is evident. The horizontal solid lines 
show the confidence levels corresponding to the probability of 
detecting a peak due to a genuine signal equal to 90% and 99%. 
The inset shows the power spectrum of the window function 
for a sine wave signal of amplitude 1 m s _1 , sampled as the 
observations. The power unit are the same as the main figure. 

segments were used in the final analysis. The internal veloc- 
ity error of a spectrum is calculated as the error of the mean 
velocity of all segments used for the analysis. The parameters 
for the instrument profile modelling were determined by using 
fast rotating featureless stars, the standard radial velocity star 
t Cet and the planet-bearing 5 1 Peg a nd p CrB stars. Detailed 
description of the technique is given in Deside ra et al. 1 0003). 
Fig. 2 shows, as an example, the radial velocity time series ob- 
tained on 2001/01/06. In this case, the radial velocity internal 
error and the r.m.s. of the data set are 1 .48 m s and 3.5 m s , 
respectively. A 21-min pulsation pattern is clearly visible. 

We have u sed the modified Scargle-Lomb periodogram 
(IScargle II 19821) in order to perform the spectral analysis. The 
result obtained by combining all the six observing nights is 
shown in Fig. 3, where an excess of power around 1 mHz is 
clearly visible. The time scale gives a spectral resolution of 
1.62 /iHz. The calibration of the Scargle's power in terms of 



lable 2. Prominent peaks present in the penodogram or 
Procyon A detected by SARG with a confidence level larger 
than 90%: bold characters indicate the frequencies consistent, 
within our spectral resolution of 1.62 //Hz, with those found by 
Martic et al. I d2004l) with a spectral resolution of 1 .79 //Hz. 

Frequencies (//Hz) 
230.2 233.7 238.4 245.4 596.5 776.7 822.8 



amplitude squared was performed in such a way that a syn- 
thetic sine wave of 1 m s _1 , sampled as the observations, would 
yield in the periodogram a cluster of peaks, the largest of which 
with a power equal to 1 m 2 s~ 2 . We then computed the statis- 
tical significance of the detected peaks in terms o f the "false 
alarm probability" F, as defined by Scarele ( 1982), and adopt- 
ing the prescription of Home & Baliunas ( 1986). F is the prob- 
ability that a spike (i.e. power in one frequency bin) is due to 
noise. Consequently the quantity (1-F) is the probability of de- 
tecting a peak due to a genuine signal. The power levels (confi- 
dence levels) concorresponding to the value (1-F) equal to 90% 
and 99%, equivalent to a false alarm probability of 10% and 
1%, respectively, are reported in Fig. 3. In Table 2 we listed 
all the individual frequencies detected with a confidence level 
grea ter than 90%; the f requencies consistent, with those found 
by Marti cet al. I (12004) are shown in bold-face. The additional 
frequencies of 302.0 ti Hz, 428.0 uHz and 8 08.8 //Hz are also 
consistent with those of Martic et al. (2004), but have been de- 
tected with a confidence level smaller than 90%. By compar- 
ing the p ower spectrum show n in Fig. 3 with the results pre- 
sented in Ma rtic et al. I (119991) . it is worth noting that most of 
the strongest peaks in our periodogram hav e power of about 0.9 
m 2 s~ 2 , a value greater than that found by Mart ic et al. I (fl999) 
in the analysis of the observations carried out on November 
1997 using the ELODIE spectrograph, which ranges from 0.8 
m 2 s~ 2 in the periodograms obtained for two subsets of data 
to 0.5 m 2 s -2 in the periodogram of the entire data set. On the 
other hand the high frequency mean white noise level in our 
periodogram, evaluated in the frequency r ange 1.8-3 mHz, is 
0.049 m 2 s -2 , while the values reported bv lMartic et al. 1 (11999) 
range from 0.014 to 0.030 m 2 s -2 , depending on the differ- 
ent set of data a nalyzed. The power of most of the strongest 
peaks reported by Martic et al. T d2004h referring to the ELODIE 
and the combined ELODIE/ AFOE data sequences colleced on 
January-Febraury 1999 is about 0.16 - 0.20 m 2 s~ 2 , while that 
reported for the AFOE data is about 0.5 - 0.6 m 2 s -2 ; the mean 
noise levels are quite different in the three sets of data, ranging 
from 0.0256 to 0.0064 m 2 s -2 . 

In order to estimate the large frequency separation Avo in 
the region of the excess of power detected in the periodogram, 
we CLEAN-ed the spectrum by the window function (fo r the 
details of the adopetd procedure see Rob erts et al. I I 1987) and 
successively we applied to th e CLEAN-ed power sp ectrum the 
comb-response (CR) method jKieldsen et al. 119951) which is a 
generalization of the power spectrum of a power spectrum and 
consequently allows us to search for any regularity in a spectral 
pattern. In the case of solar-like oscillations it allows us to de- 
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Fig. 4. Comparison between modified comb responses for 
the two centrai frequencies 738 //Hz and 1321 //Hz of the 
power spectrum of Procyon A data. The responses are 
scaled to the maximum peak obtained for the frequency 

v = 1321 //Hz at Av = 47.8 //Hz. 



(MHz) 



Fig. 5. The first-order spacings over several frequency 
ranges, as determined from comb analysis. The line shows 
the mean value of 56 //Hz. 

tect their typical comb-like structure. In particular, a peak in the 
CR at a particular spacing Avo indicates the presence of a reg- 
ular series of peaks in the power spectrum, centered at the cen- 
tral frequency vo and having a spacin g of Avn/2. We ado pted 
the modified comb-response function ( Marti c et al. l2004l) : 



C(v , Av ) = IU Q [S(vo + & + 6)S(v - if* + 6W> 



■AV 



(1) 



where 6 = ±(i mod 2)Dq according to mode degree I — or 
I = 1 at the central frequency vo and a-, = 1 for i — 1 or 2, 
a, = 0.5 for i = 3 or 4. For each central frequency vo, al- 
ternatively considered as / = and / = 1, we searched for 
the maximum CR in the intervals 20 < Avo < 80 fiHz and 
0.3 < Dq < 2. Fig. |4] shows examples of the modified CR 
computed for the two central frequencies 738 //Hz and 1321 
//Hz. In Fig. [5] we show the variation of Avo in the frequency 
range 500 < vo < 1500 //Hz, as determined from CR analy- 
sis. The arithmetic mean of this spacing is 56 + 2 //Hz, a value 
of the large separ ation in good agre ement with both theoreti- 
cal estimates dChabover etal.ll999l) and previous observations 



( Martic et al . 2004 ) . 1 he error on A vo has been estimated as the 
standard deviation of the Avo's computed from CR analysis. 

The determination of Avo is not strongly sensitive to the 
choice of Do. 

4. Conclusions 

The analysis of the periodogram of Procyon A obtained by 
SARG shows a p-mode spectrum characterized by a large 
freq uency separation A vn — 56 + 2i/Hz in ag reement with 
Cha bover et al I i 19991) andlMartic et al. I (|2004). Contrarily to 
what found bv lMatthews etal. M2004I) our results strongl y sup- 
port t he idea that the excess of power detected by M artic et al. I 
(2004) in the range from 0.5 up to 1 .5 mHz is caused by p-mode 
oscillations. In particular we think that the SARG spectrograph 
can efficiently be used in future multi-site observing campaigns 
for asteroseismology. 
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